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ABSTRACT: Beam theory analysis was applied to predict delamination growth in DCB 
specimens reinforced in the thickness direction with pultruded pins, known as Z-fibers. 
The specimen arms were modeled as cantilever beams supported by discrete springs, 
which were included to represent the pins. A bi-linear, irreversible damage law was used 
to represent Z-fiber damage, the parameters of which were obtained from previous 
experiments. Closed-form solutions were developed for specimen compliance and 
displacements corresponding to Z-fiber row locations. A solution strategy was 
fonnulated to predict delamination growth, in which the parent laminate mode I fracture 
toughness was used as the criterion for delamination growth. The solution procedure was 
coded into FORTRAN 90, giving a dedicated software tool for performing the 
delamination prediction. Comparison of analysis results with previous analysis and 
experiment showed good agreement, yielding an initial verification for the analytical 
procedure. 


KEYWORDS: Delamination prediction, DCB specimen, strain energy release rate, 
Simple beam theory, Z-fibers. 
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LIST OF SYMBOLS 


a; Distance from spring at location i to beam root. 
a 0 Length of non-reinforced laminate. 

a z Z-fiber spacing. 

A p Z-fiber cross sectional area, 

b DCB Specimen width. 

C DCB Specimen compliance, 

d Damage parameter. 

E p Z-fiber axial tensile modulus. 

E Flexural modulus of DCB specimen arm with no Z-fiber reinforcement. 

E zp Flexural modulus of DCB specimen arm with Z-fiber reinforcement. 

F c Critical load for onset of Z-fiber pull-out. 

F z Beam forces in z-axis. 

G Strain energy release rate. 

Gi c Critical mode-I strain energy release rate, 

h DCB specimen ann height. 

I Moment of inertia of DCB specimen arm with no Z-fiber reinforcement. 

I zp Moment of inertia of DCB specimen arm with Z-fiber reinforcement, 

k; Current Spring stiffness. 

k 0 Original spring stiffness. 

L Total delamination length. 

L p Z-fiber length. 

M Bending moment, 

n Mode-I compliance calibration exponent. 

P Applied load. 

P c Critical load for onset of delamination growth in a DCB specimen, 
r Number of Z-fibers per row. 

R Spring reaction load. 

U Stored elastic strain energy. 

Zf Z-fiber row displacement at complete pull-out. 

z; Current Z-fiber row displacement. 

z 0 Z-fiber row displacement at onset of pull-out. 

6 DCB specimen load-point displacement. 

A Mode-I crack length extension factor. 

Adev Mode-I crack length extension factor during developing phase of crack growth. 

A ss Mode-I crack length extension factor during steady-state phase of crack growth. 

o c Critical stress for onset of Z-fiber pull-out. 
a) Ratio of Z-fiber spacing to Z-fiber insertion depth. 
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INTRODUCTION 


A method for enhancing fracture toughness of laminated composite structures uses pins 
inserted normal to the laminate plane to increase delamination resistance. The pins, 
known as Z-fibers, are inserted into the prepreg laminate using an ultrasonic hammer 
prior to the curing process, resulting in a field of pins embedded normal to the laminate 
plane as illustrated in Fig. 1 . A cross section of laminate containing Z-fiber reinforcement 
is pictured in Fig. 2. It can be seen that the reinforcing fibers are not exactly normal to the 
laminate plane. This is a consequence of the fibers rotating during the insertion process. 

Z-fibers are typically 0.28-mm to 0.5-mm in diameter. Standard areal densities 
range from 0.5% to 4%. The Z-fibers are provided by the manufacturer, Aztex®, in a 
low-density foam preform, which acts to stabilize orientation of the pins during the 
insertion process [1-3]. Typical pin materials include boron and carbon fibers embedded 
in a polymer matrix. 

Z-fibers act to enhance mode I delamination resistance by providing a direct 
closure traction to the opening of delaminating surfaces. A study of delamination growth 
in IMS/924 double cantilever beam (DCB) specimens [4], reinforced with carbon/BMI Z- 
fibers, demonstrated that a significant increase in Mode I fracture toughness is achieved 
as a result of the through thickness reinforcement [5]. It was found that an areal density 
of 0.5% and 2% Z-fibers resulted in a four-fold and seventeen-fold increase in the non- 
reinforced fracture toughness of the laminate respectively. 

Behavior of delamination growth in DCB specimens containing Z-fibers is 
defined as a three-phase process in the current work, although others consider it a two- 
stage process [6-7]. Firstly, the delamination front reaches the Z-fiber field (Stage 1, Fig. 
3). Delamination continues through the field of reinforcement, causing Z-fibers to bridge 
the delamination. Initially, the number of bridging Z-fibers increases as the delamination 
grows, resulting in an increase in critical fracture load. This is termed the developing 
growth phase and is labeled Stage 2 in Fig. 3. Finally, delamination growth reaches a 
steady-state, during which a new bridging Z-fiber row is created as another is completly 
pulled out from the material. The net effect is a constant number of bridging Z-fiber 
rows as delamination growth continues. This is termed Steady-state growth (Stage 3, 
Fig-3). 

A number of researchers have performed numerical and analytical studies of 
delamination growth in DCB specimens, reinforced either by stitching or Z-fibers. Jain 
and Mai [8] developed an analytical model of a stitched DCB specimen, in which the 
bridging stress provided by the stitches was represented as a distributed load over the 
bridged zone of the specimen arm. The analyses were used to predict the increase in 
delamination resistance associated with the through-thickness reinforcement. 

During an investigation into the fracture behavior of debonding in integrally 
stiffened panels, Glaessgen et al [9] performed finite element analyses of stitched DCB 
specimens. In this work, plane stress plate elements were used to represent the skin and 
stiffener materials, while the stitches were modeled as discrete, nonlinear, fastener 
elements, the stiffness of which was obtained from experiment. The virtual crack closure 
technique (VCCT) [10] was used to compute energy release rate values following each 
analysis. The analysis was used to study the effect of stitching parameters on energy 
release rate of the reinforced specimen. 
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Cartie and Partridge [11] performed finite element analyses of unidirectional, 
carbon fiber/epoxy DCB specimens containing Z-fibers. Plane strain elements were used 
to represent the laminate material. Single Z-fiber pull-out tests were performed prior to 
the analyses, to determine the bridging law exhibited by the Z-fibers as they are loaded 
under mode I conditions. Loading was applied at nodal locations corresponding to the 
position of each Z-fiber row, in order to represent the bridging action exhibited by the 
through thickness reinforcement. Load-displacement plots and delamination resistance 
curves (R-curves) were generated following the analyses. 

In ref. 12, Robinson and Das performed DCB tests on IM7/8552 tape laminates 
with a unidirectional stacking sequence and reinforced with Z-fibers. An analysis similar 
to ref. 8 was developed to predict delamination growth in the DCB specimens. The 
analysis was conducted in two parts. First, the developing phase of delamination growth 
was analyzed. In the second stage of the analysis, delamination growth corresponding to 
the steady-state phase was predicted. The bridging stresses exhibited by the bridging Z- 
fibers were modeled as a distributed load along the bridged section of the DCB specimen 
arms. Following the analysis, compliance calibration and modified beam theory 
techniques [4] were used to calculate apparent fracture toughness values. 

The aim of the current work was to develop an analysis technique capable of 
predicting delamination growth behavior of DCB specimens reinforced with Z-fibers. 
The resulting technique modeled specimen arms as cantilever beams with discrete springs 
used to represent Z-fiber rows. Closed-form solutions were derived for specimen 
compliance and Z-fiber displacement. The solution procedure was coded into 
FORTRAN 90, enabling an efficient execution of the analysis. The analysis can be used 
to study the effect of Z-fiber insert length, areal density and diameter on the apparent 
fracture toughness of the reinforced specimen. Furthermore, the analysis is coded to 
enable the use of alternative bridging laws for representing Z-fiber failure. Because Z- 
fibers are represented using springs, the unstable nature of delamination growth is 
captured by the analysis, unlike the analysis of refs. 8 12, in which the bridging stresses of 
the Z-fibers are represented as a distributed load. The analysis technique may also be 
used to design Z-fiber reinforced DCB specimens. Preparation for DCB tests can be 
eased by using the analysis to estimate the maximum load and displacement that will be 
experienced during a test. 


Analysis Procedure for Predicting Delamination Growth in Z-Fiber Reinforced 
DCB Specimens 

The analysis procedure used for predicting delamination growth in Z-fiber reinforced 
DCB specimens is comprised of several main steps. First, the specimen is modeled as 
cantilever beams with springs representing Z-fiber rows. Closed-form solutions for the 
specimen compliance and displacements are then derived. A damage law is applied to 
the springs, in order to represent the behavior of Z-fibers as they are removed from the 
parent laminate during a DCB test. Linear elastic fracture mechanics is used to calculate 
the strain energy release rate, G, in the specimen for a given load-point displacement. 
The value of G is compared with the mode-I critical strain energy release rate, Gi c , of the 
parent laminate material, and delamination growth is assumed to take place when G is 
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greater than or equal to Gi c . Finally, a solution strategy was developed to allow 
delamination growth prediction in a Z-fiber reinforced DCB specimen. The remainder of 
this section details the five main steps of the analysis. A description of the data reduction 
schemes that are used to calculate Gi c from the analysis data is given at the end of the 
section. 


Z-Fiber Reinforced DCB Specimen Model 

Z-fiber reinforced DCB specimens are modeled using cantilever beams to represent the 
specimen arms. Discrete springs are used for representing the Z-fibers that may bridge 
the delamination. Figure 4 illustrates how the 2-dimensional model represents an actual 
DCB specimen, where three rows of Z-fibers are shown to be bridging the delamination. 

Displacement Equations of a Z-Fiber Reinforced DCB Specimen 

Closed-form solutions are derived for the displacement of DCB specimens with 
delaminations bridged by one, two and three Z-fiber rows. Finally, the solutions are 
generalized to accommodate for any number of bridging Z-fiber rows. 

It is assumed that the upper and lower DCB specimen arms are identical and 
displacement of each arm is symmetrical about the neutral axis of the specimen, as 
illustrated in Fig. 4. Consequently, a single specimen arm is modeled throughout the 
analysis. The displacement equations are doubled, yielding an estimate of the total DCB 
specimen displacements. 

A schematic of a specimen arm model with one bridging Z-fiber row is given in 
Fig. 5a. The free body diagram of the arm is given in Fig. 5b. Z-fiber rows were 
represented by discrete springs positioned normal to the x-y plane. The length of the 
spring is equivalent to the Z-fiber insertion depth in one arm. Studies have shown that in- 
plane elastic properties of a laminate may suffer some degradation as a result of the Z- 
fiber insertion process [13-15], Consequently, the flexural modulus of the laminate 
containing Z-fibers was treated separately from that of the remaining laminate. Using the 
free body diagram in Fig. 5b, the equilibrium equations of the arm are: 


2F,-0-P-R„-^-(2k 1 ) .-. 

2M-0-PL-M o -|(2k,)a 1 


R.-P-yPD 

M 0 -PL-^(2k,) a 


( 1 ) 


In the above and remaining equilibrium equations, the spring displacements, z, are halved 
and the stiffness values, k, are doubled to accommodate for the fact that half the total Z- 
fiber is modeled. 

The reaction load, Ri of a spring at location ‘i’ is assumed to be in direct 
proportion to the displacement of the spring in the z-axis, z h by the spring stiffness, k, 
such that: 


R; = 



( 2 ) 
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z 

Following simple beam theory analysis, the spring displacement,-^-, and load-point 
displacement, — of one half of the DCB specimen are expressed as: 


z, Paj 2 (3L-aj) 
2 = 6E zp I zp + 2k ia ? 


2 

Where E and E zp are the flexural modulus of the non-reinforced and Z-fiber reinforced 
sections of the arm respectively. I is the moment of inertia of the non-reinforced 
specimen arm and I zp is the moment of inertia of the specimen arm containing Z- fibers. 
These moment of inertia values are equal in all cases studied in the current work. L is the 
arm length, measured from the loading point to the crack tip. ai is the distance between 
the Z-fiber and crack tip. The parameter a 0 is the length of non-reinforced laminate, 
where the region of Z-fiber reinforcement is assumed to begin one Z-fiber spacing away 
from the first Z-fiber row. P is the load applied to the specimen, kj is the spring stiffness 
Assuming the upper and lower arms of the specimen are identical and the 
delamination is located at the specimen mid-plane, Eqns.3 and 4 are doubled to yield the 
Z-fiber row displacement and load-point displacement of the complete DCB specimen: 


L 3 -a 3 


3EI 


3E I 

zpzp 


kjZja^aj - 3L) 
6E I 

zp zp 


( 4 ) 


Pa 2 (3L - a,) 

3EX+VF 


6, = 2P 


L 3 -a 3 


3EI 


3E I 

zp zp 


k^a^a, - 3L) 
3E I 

zp zp 


( 5 ) 

( 6 ) 


The subscript on the left-hand side of Eqn.5 indicates displacement of the specimen at a 
distance, ai, from the arm root, as illustrated in Fig. 5a. The subscript on the left-hand 
side of Eqn.6 represents the number of Z-fibers bridging the specimen. 

Specimen compliance, C, is calculated by dividing the load-point displacement by 
the applied load, P. Therefore, compliance of the DCB specimen with a single bridging 
Z-fiber is found from Eqn.6: 


C, = — = 2 
1 P 


L 3 - a 3 


3EI 


3E I 

zpzp 


kjZja^a, - 3L) 
3PE I 

zp zp 


( 7 ) 


Again, the subscript on the left-hand side of Eqn.7 refers to the number of bridging Z- 
fiber rows. 
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The analysis is now extended to calculate the displacements of a DCB specimen, in 
which the delamination is bridged by two Z-fiber rows. Figure 6a illustrates the 
specimen arm. A free body diagram of the arm is also shown in the figure. The 
equilibrium equations of the arm are: 


2F,-0-P-R Q -|(2k,)-| 

2M-0-PL-M„-|(2k,)a, 


(2k.) R 0 -P-(|(2k,)-|(2k 2 )j 

-y(2k> 2 .-. M„-PL-(|(2k l )a l - 



( 8 ) 


Beam theory analysis leads to solutions for the Z-fiber row displacements*, zi and Z 2 , at 
locations 1 and 2 respectively along the beam. 


Paf(3L - a,) + z 2 k 2 a 1 ’(a 1 - 3a 2 ) 
3E zp I zp + k iai 3 

Pa 2 (3L - a 2 ) + - 3a 2 ) 

3E zp I zp +k 2 a 2 


(9) 

( 10 ) 


Compliance of a DCB specimen with two bridging Z-fiber rows is: 


^-2 

3 y 3 3 ' 

a o , L “ a o 

_l_ 

k i z i a i ( a i 3L) + k 2 z 2 a 2 (a 2 3L) 

P 

3EI 3E zp I zp 

T 

3PE I 

zp zp 


( 11 ) 


The schematic of a specimen arm bridged by three Z-fiber rows is given in Fig. 6b. The 
free body diagram of the arm is also shown and the corresponding equilibrium equations 
are: 


2 F,-0-P-R o -|(2k,)-|(2k 2 ) -|(2k 3 ) 



Z-fiber row displacements, zi, z 2 and z 3 , at locations 1, 2 and 3 respectively along the 
beam are: 


The remaining solutions are doubled to reflect displacements and compliance of the entire DCB 
specimen. 
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(13) 


z i = 


Pa, (3L - a,) + z 2 k 2 af(a l - 3a 2 ) + z 3 k 3 af(a 1 - 3a 3 ) 


3E zp I zp + k,a. 


Pa 2 (3L - a 2 ) + z^a^a, - 3a 2 ) + z 3 k 3 a 2 (a 2 - 3a 3 ) 
Z ' = 3E zp I zp + k 2 a 2 

Pa 3 (3L - a 3 ) + z,k 1 a 1 2 (a 1 - 3a 3 ) + z 2 k 2 a 2 (a 2 - 3a 3 ) 
23 = 3E zp I zp + k 3 a 3 

Compliance of the DCB specimen with three bridging Z-fiber rows is: 


^-2 

3 y 3 3 " 

a o , L - a o 

_L 

kjZja^a, - 3L) + k 2 z,a 2 (a 2 - 3L) + k 3 z 3 a 3 (a 3 - 3L) 

P 

.3EI 3E zp I zp 

T 

3PE I 

zpzp 


(14) 

(15) 


(16) 


The analysis is extended to included any number of Z-fiber rows bridging the 
delamination of the DCB specimen. Figure 6c shows one arm of a DCB specimen 
bridged by an arbitrary number of Z-fiber rows, n. The free body diagram of the arm is 
also given in Fig. 6c. Following the diagram, the equilibrium equations for the arm are: 


2e z = 0 = p-r o - 




(17) 


2M = 0 = PL-M o -2|(2k i )a i M o = PL - J|(2k i )a i 

i= 1 i= 1 

Displacements in the z direction, corresponding to each Z-fiber row location in the DCB 
specimen are: 
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z l = 


3E zp I zp + kjSj 


Pa 2 (3L - a,) + z 2 k 2 a 2 (a, - 3a 2 ) + z 3 k 3 a l ’(a 1 - 3a 3 ) + • • • 

+ Z n-l k n-i a i 2 ( a i - 3a n-l) + Z „ k A( a i " 3a „) 


Z 2 = 


3E zp I zp + k 2 a 2 


Pa 2 (3L - a 2 ) + z^a 2 ^ - 3a 2 ) + z 3 k 3 a 2 (a, - 3a 3 ) + — 

+ Z „4 k „-I a 2( a 2 - 3a „-l) + Z n k n a 2( a 2 " 3a n ) 


3E zp I zp + k 3 a 3 


Pa 3 (3L - a 3 ) + z i k 1 af(a 1 - 3a 3 ) + z 2 k 2 a 2 (a, - 3a 3 ) + — 

+ Z n-l k „-i a 3( a 3 - 3a „-l) + Z n k n a 3( a 3 " 3a n ) 


( 18 ) 


Z „-l = 


3E I +k ,a 3 


Pa n-l( 3L - a n4 ) + z i k i a f( a i - 3a n -1 ) + Z 2 k 2 a 2( a 2 - 3a n-l) 

z 3 k 3 a 3 (a 3 - 3a n4 ) + • • • • +z n k n a 2 _ 1 (a n ^ 1 - 3a n ) 


Pa n(3L -a n ) + Z ! k i a i 2 ( a i - 3a „) + Z 2 k 2 a 2( a 2 " 3a „) 

z 3 k 3 a 3 (a 3 - 3a n ) + • • • • +z n . 1 k n4 a 2 . 1 (a n . 1 - 3a n )_ 

The displacements are solved simultaneously, yielding the solutions for each Z-fiber row 
displacement. 

Compliance of the DCB specimen bridged by any number, n, Z-fiber rows is 
given by: 


z„ = 


3 E zp i zp + k n a : 


c 


n 



L 3 - a 2 


3E I 

zp zp 


3PE I 

zp zp 


2 k i z i a ?( a i- 3L ) 

. i=l 


(19) 


Z-Fiber Pull-out Bridging Law 

A bi-linear, softening constitutive model [16] was used to represent the damage incurred 
in Z-fiber rows, as they pull from the laminate during a DCB test. A typical illustration 
of the history from a Z-fiber pull-out test is given in Fig. 7. Upon initial loading, the Z- 
fiber row exhibits a constant stiffness, until a critical displacement, z 0 , is achieved, after 
which the stiffness decreases linearly with displacement. Load displacement responses 
following pull-out tests on single Z-fibers have exhibited similar damage behavior [11, 
12 ]. 

Using this damage model, the stiffness of the Z-fibers is calculated according to the 
following conditions (Fig. 7): 


k o <= Z , £Z o 

k; = \ (l-d)k Q <=z Q <z ; <z f 
0 <= z. > z f 


( 20 ) 
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Where z 0 is the Z-fiber row displacement corresponding to the onset of damage, Zf is the 
displacement at which the Z-fiber row stiffness becomes zero, corresponding to complete 
pull out. It is usual for Z-fibers to be removed from one side of the DCB specimen only. 
Therefore, Zf is assumed to be half the total Z-fiber insert length. The damage parameter, 
d, is given by [16]: 


d _ Z f( Z i~ Z o) 
Zi(Zf-Zo) 


( 21 ) 


The initial Z-fiber row stiffness, k 0 , is calculated using the following relation: 


k = r 


E p A p 


( 22 ) 


Where Ei, A p and L p are the axial tensile modulus, cross sectional area and length of a 
single Z-fiber respectively, r is the number of Z-fibers per row. 

The Z-fiber row displacement corresponding to the onset of damage, z 0 , is 
calculated from the row stiffness, k 0 , and the peak load, F c , measured from a single Z- 
fiber pull out test, (Fig. 7) such that: 


z 


o 



(23) 


The Z-fiber displacement corresponding to complete removal from the specimen, Zf, is 
equal to the insertion depth into one of the DCB specimen arms. 

Strain Energy Release Rate Calculation 

Strain energy release rate values were calculated by repeating an analysis of a DCB 
specimen containing delaminations of neighboring lengths, as illustrated in Fig. 8a. The 
same load-point displacement, 6, is applied in both analyses. The strain energy, Ui, 
stored in the specimen with delamination of length Li is equal to the area OAC under the 
load-displacement plot illustrated in Fig. 8b. The stored elastic strain energy, U 2 , in the 
specimen with the neighboring delamination of length, L 2 is given by the area OBC under 
the load-displacement plot in Fig. 8b. Therefore, the strain energy release rate, G, for 
delamination growth from Li to L 2 is given as: 


J_dU _ 1. U.-U, 
b dL b L 2 -L 1 


(24) 


Where b is the DCB specimen width. 

The strain energy release rate for a delamination length of Li is estimated from 
Eqn.24. This calculation of G becomes more accurate as the difference between 
neighboring delamination lengths, Li and L 2 , decreases. Consequently, a convergence 
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study was undertaken to determine the value of (L 2 -Li) required to yield accurate strain 
energy release rate estimates. This study is detailed in the model verification section. 

Delamination Growth Solution Strategy 

A strategy was formulated to apply the analysis detailed above, enabling prediction of 
delamination growth through a Z-fiber reinforced DCB specimen. 

Firstly, the specimen configuration is prescribed, including specimen dimensions, 
elastic properties, areal density of Z-fibers, number of Z-fiber rows and Z-fiber insertion 
length. An initial load-point displacement, 6, is applied and the corresponding specimen 
compliance is calculated. The following relation is then used to determine the applied 
load, P: 


P = 


h_ 

C 


(25) 


P, is then used to calculate the Z-fiber row displacements. If more than one row bridges 
the delamination, Gaussian elimination with partial pivoting is used to solve the 
simultaneous equations. The method described previously is then used to calculate the 
strain energy release rate corresponding to the displacement increment. This value is 
then compared with the critical mode I strain energy release rate, Gi c , of the parent 
laminate. Delamination is assumed to occur when G = Gi c and thus a new delamination 
length increment is chosen at this stage. If G < Gi c , the displacement increment is 
repeated using a larger value of load-point displacement. Once a Z-fiber row 
displacement reaches the critical value for damage onset, the bi-linear bridging law is 
implemented to adjust the Z-fiber row stiffness. During this stage of the solution 
procedure, each displacement increment is repeated with the updated Z-fiber row 
stiffness values until no further change in stiffness is incurred within that displacement 
increment. Once equilibrium of the row stiffness values has been achieved, the energy 
release rate calculations are performed and the remainder of the analysis continues. A 
flow chart summarizing the solution procedure, as implemented into the FORTRAN 90 
code, is presented in Fig. 9. 


Data Reduction Schemes 

The load-point displacement, applied load and corresponding delamination lengths 
calculated from the analysis are reduced in a similar manner to data obtained from an 
actual DCB test. Compliance calibration and modified beam theory data reduction 
schemes [4] were selected for calculating an effective fracture toughness of the Z-fiber 
reinforced DCB specimen. The corresponding expressions for fracture toughness using 
the compliance calibration and modified beam theory methods are given in Eqn. 26. 


G IC = 


G IC = 


nPA 

2bL 

3PA 

2b(L +A) 


(Compliance calibration) 
(Modified beam theory) 


(26) 
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Where P c and 5 C are the critical load and displacement respectively for the onset of 
delamination growth. 

These data reduction schemes require that the term, n, is determined as the slope 
from a plot of the log of compliance versus log of load-point displacement. A crack 
length extension factor, A, is determined as the intercept from a plot of the cube root of 
compliance versus delamination length. 


Model Verification 

Verification of the current modeling approach was undertaken in four stages. Firstly, a 
convergence study was conducted on the method for calculating strain energy release rate 
values using Eqn.24. The current analysis was used to model two independent specimen 
configurations available in the literature. Finally, a series of analyses were conducted to 
demonstrate that the current analysis captures unstable delamination growth, when the 
spacing of Z-fiber rows is sufficiently large. 

Convergence Study on Non-reinforced DCB Specimen [17] 

A convergence study was conducted to determine the difference in neighboring 
delamination lengths, dL, required to yield accurate strain energy release rate values, as 
calculated using Eqn.24. The DCB specimen discussed in reference 17 was used in the 
convergence study. The specimen contained no Z-fiber reinforcement. Other specimen 
details are presented in Table 1. Critical load and load-point displacement values were 
estimated for delamination lengths of 50-mm and 99-mm. The analyses were repeated 
using different values of dL. The estimated critical load and displacement values for 
delamination growth were normalized by the corresponding experimental values and 
plotted as a function of dL/delamination length, L. The results indicate that neighboring 
analyses differ by less than one percent once the difference in neighboring delamination 
lengths is equal to or less than 0.2% of the delamination length, Li The results for the 50- 
mm and 99-mm cases are presented in Figs. 10a and 10b respectively. In both cases 
convergence is assumed when dL/L < 2xl0" 3 , and therefore this value of dL was used for 
the difference in neighboring crack lengths. To ensure convergence was achieved during 
analysis of Z-fiber reinforced DCB specimens, the difference in neighboring 
delamination length was set to 0.1% of the Z-fiber spacing. 

Analysis of Z-Fiber Reinforced IM7/8552 DCB Specimens [12] 

Robinson and Das conducted DCB tests on specimens manufactured from IM7/8552 tape 
laminate with a stacking sequence of [0] 48 , yielding a specimen arm thickness of 3 -mm. 
A field of carbon/BMI Z-fibers was inserted into the DCB specimens as illustrated in 
Fig. 1 1 . Z-fibers were placed at 2 -mm intervals, giving an areal density of 1.5%. An 
analysis method was used to predict delamination growth in the DCB specimens. Critical 
strain energy release rates were calculated using the data reduction schemes described 
earlier. 

The load displacement responses predicted following the original and current 
analyses are presented in Fig. 12 and show reasonable agreement. Data reduction from 
the current analysis was carried out in order to calculate an apparent toughness of the 
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specimen. Compliance calibration and modified beam theory plots from the current 
analysis data are presented in Figs. 13 and 14 respectively. Both curves exhibit linear 
regions that correspond to the developing and steady-state phases of delamination 
growth, as indicated in the figures. The steady-state phase of delamination growth 
yielded the highest value of the exponent, n, and the lowest crack length correction 
factor, A, in comparison to the developing stage. Similar observations were made by 
Robinson and Das [12]. Compliance calibration and modified beam theory were used to 
calculate the apparent fracture toughness and were found to produce similar toughness 
values. A plot of the R-curve as predicted from the current analysis is presented in 
Fig. 15. Also in the figure is the R-curve prediction from the original analysis [12]. The 
toughness values shown from the current analysis were calculated using the compliance 
calibration method, and a value of n = 3.19. 

Predictions from the two analyses agree well. However, the data reduction 
schemes, which were developed for non-reinforced laminates may produce invalid global 
toughness values because the damage zone is significant relative to the specimen 
dimensions (as the damage zone spans the bridging Z-fibers being removed from the 
specimen) and in the developing stage, delamination growth is not self-similar. Because 
self-similar delamination occurs only during the steady-state phase of delamination 
growth (when the number of bridging Z-fibers is constant), only the data corresponding 
to this delamination stage is recommended for use with existing data reduction schemes. 
Similar recommendations were made by Robinson and Das [12], 

Analysis of Z-Fiber Reinforced IMS/924 DCB Specimens [17] 

The current analysis technique was used to model DCB tests conducted on IMS/924 
carbon tape specimens with a unidirectional stacking sequence, reinforced with 
carbon/BMI Z-fibers [17]. A specimen configuration containing 0.5% areal density and 
one containing 2% areal density of Z-fibers was modeled, details of which are presented 
in Tables 3 and 4 respectively. The bridging law for each case was obtained from single 
Z-fiber pull-out tests [17]. 

The load displacement response estimated from the analysis of both specimen 
types is given in Figs. 16 and 17 respectively. A comparison with the experimental load 
displacement traces is given for each case and shows a reasonable agreement between 
experiment and analysis in the 0.5% areal density case (Fig. 16). The developing section 
of the load displacement traces in the 2% areal density case calculated from the analysis 
does not agree well with the experiment. This disagreement, however, does reduce in the 
steady-state stage of delamination growth (Fig. 17). A selection of delamination lengths 
corresponding to the analytical load-displacement traces is also given. Figures 18 and 19 
present plots of the load-point displacement versus delamination length from the analysis 
and corresponding experiment for the 0.5% and 2% areal density cases respectively. The 
first case shows good agreement with experiment (Fig. 18), indicating the analysis yields 
a reasonable prediction for delamination growth. In the 2% areal density case (Fig. 19), 
the analytical displacement versus delamination length curve only agrees with experiment 
after steady-state growth has been reached. The disparity is likely due to the assumptions 
made in the analysis and the variability of the Z-fiber region in the experiment, as 
parameters such as insertion depth and angle of Z-fibers have been shown to greatly 
affect the resulting delamination response of DCB specimens [5]. Furthermore, uneven 
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delamination growth is likely to take place, which will influence the delamination length 
readings taken during the experiment. This disparity in load-point displacement versus 
delamination length curves may also explain the difference in the developing regions of 
the load displacement traces from analysis and experiment for the 2% areal density case. 

Again, the compliance calibration and modified beam theory techniques were 
used to reduce the analytical data. Consider first the compliance calibration plots from 
analyses of both specimens shown in Fig. 20. In accordance with the previous example, 
the steady-state phase of growth yields the highest values of n. The corresponding 
modified beam theory curves for both specimens are given in Fig.21. Again, the values 
of crack length extension factor are highest during the developing stage of delamination 
growth. A summary of the findings from the compliance calibration and modified beam 
theory plots are given in Figs. 20 and 21. 

Because the developing phase of delamination growth is not self-similar, the 
apparent fracture toughness was calculated for the steady-state phase only. The 
experimental data corresponding to the steady-state phase of delamination was used in 
the current work with the compliance calibration technique for calculating the 
experimental values of Gi c . Figure 22 presents the experimental compliance calibration 
curves for the 0.5% and 2% areal density cases. A summary of the resulting values of the 
exponent, n, is also given in the figure and shows the same trends as determined through 
the analysis. The values corresponding to steady-state delamination were used in the data 
reduction for calculating the experimental values of Gi c . The R-curves from the analyses 
for the 0.5% and 2% areal density cases are given in Figs. 23 and 24 respectively. There 
is reasonable agreement between the analytical and experimental values. The average 
experimental steady-state value of Gi c for the 0.5% areal density specimen was 1457- 
N/m. The corresponding analytical value was 1390-N/m, corresponding to a 5% error 
between experiment and analysis. The average experimental steady-state value of Gi c for 
the 2% areal density specimen was 5847-N/m and the corresponding analytical value was 
5342-N/m. The error between experiment and analysis for this case was 10%. 

The Effect of Z-Fiber Row Spacing on Delamination Growth Behavior 
If the spacing of Z-fiber rows in the direction of delamination growth is small, the direct 
closure traction provided by the bridging fibers will be similar to a distributed load across 
the bridged region of the DCB specimen. Consequently, steady-state delamination 
growth will be stable with respect to load-point displacement. However, as Z-fiber row 
spacing increases, the bridging stresses will become discrete and the resulting steady- 
state delamination will be unstable. To test the ability of the current analysis to capture 
this change in delamination behavior with Z-fiber row spacing, a DCB specimen analysis 
was run for the specimen detailed in Fig. 25a, in which the Z-fiber spacing was 1.5-mm 
(spacing to Z-fiber insertion depth ratio, (0=1.5/2=0.75). The analysis was repeated for 
the same specimen but using a Z-fiber spacing of 5-mm ((0=2.5). In order to provide a 
direct comparison of delamination behavior between the two cases, the peak stress for the 
onset of Z-fiber pull-out was kept constant. This was achieved in the following manner: 

The assumed bridging law for the case of a Z-fiber row spacing of 1.5-mm is 
summarized in Fig. 25b. The critical stress corresponding to the onset of Z-fiber pull-out, 
a c , is given by the critical load divided by the area of a unit-cell containing one Z-fiber 
[12] (Fig.25a): 
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a 


(27) 



Thus, for the 1.5-mm Z-fiber spacing example, the critical stress for the onset of Z-fiber 
pull-out was 3.4N/( 1.5mm) 2 =1.5 1 N/mnr. 

The critical stress for the 5 -mm Z-fiber spacing example was kept equal to 
1.51N/mm 2 by factoring the original critical load, P c , by the ratio of the area of the 5-mm 
square unit-cell and the 1.5-mm square unit cell, such that: 

5 2 

F = 3.4 — r- = 37.78N (28) 

c 1.5 2 

The resulting Z-fiber pull-out bridging law for the 5-mm spacing example is given in 
Fig. 25c. 

Load displacement plots from analyses of the 1.5-mm and 5-mm Z-fiber spacing 
models are given in Fig. 26. The values of load and displacement from both cases are 
expected to be the same because the critical load for the onset of Z-fiber pull-out was 
adjusted to keep the critical stress values the same. As expected, the load displacement 
response corresponding to the 5-mm spacing consists of sudden load drops. This is in 
contrast to the stable response from the example with the smaller spacing. The 
corresponding load-point displacement versus delamination length curves are presented 
in Fig. 27. In the 5-mm spacing example, steady-state delamination growth occurs via 
discrete increases for little change in displacement. This confirms that the sudden load 
drops in the load displacement response do correspond to unstable delamination growth. 

This study has shown that the current analysis method is able to capture the 
discrete delamination growth behavior of Z-fiber reinforced DCB specimens. The study 
also suggests that the representation of bridging stresses by a distributed load becomes 
less valid when the Z-fiber spacing is significantly larger than the insertion depth of the 
reinforcement. In this example significant unstable delamination growth was detected 
when the Z-fiber spacing was 250% of the insert depth of the Z-fibers. 


Discussion 

An initial verification of the current analysis for predicting delamination growth in Z- 
fiber reinforced DCB specimens has been conducted. Comparison of data from current 
analysis and a previous analysis from ref. 12 shows good agreement. This study indicated 
that the use of current data reduction schemes intended for non-reinforced DCB 
specimens is invalid when the specimens contain Z-fiber reinforcement. This is because 
delamination growth during the developing stage is not self-similar, as the number of 
bridging Z-fibers increases with delamination length. The steady-state phase of 
delamination, however, is self-similar, because the number of bridging Z-fibers is 
constant with delamination length. Subsequently, data from this stage of delamination 
growth is only recommended for use with current data reduction strategies. Similar 
conclusions were drawn in the study of ref. 12. This strategy for data reduction was 
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applied to experimental data from ref. 17 and corresponding analyses were undertaken. 
The resulting comparison between the current analysis and experiment is acceptable, 
although there is some scatter in the experimental compliance calibration curves. 
However, the experimental data exhibited the same trends as the analysis results, in that 
the exponent, n, corresponding to steady-state delaminating growth was larger than the 
developing stage value. 

The Z-fiber spacing study demonstrated that the current analysis is able to capture 
unstable delamination growth, which in this instance occurs when the spacing of Z-fiber 
rows is greater than the insertion depth of the Z-fibers. 

The studies undertaken using the current analysis have provided an initial 
verification of the accuracy of the analytical method. The analysis, however, does suffer 
a small number of limitations, such as the inability for representing a Z-fiber insertion 
angle different from 90 degrees to the laminate plane. The analysis does not 
accommodate for displacement of the DCB specimen arms due to shear, which will 
become significant in specimens with thick arms. The analysis is also limited by the 
assumption made in the beam theory that the beam displacements are small relative to the 
span. The accuracy of the analysis is also limited by the accuracy of the empirical 
damage law used to represent Z-fiber failure. In summary, the current analysis shows 
good agreement experimental results but further study will be needed to overcome 
current limitations. 


Concluding Remarks 

Simple beam theory was used to derive closed-form solutions for the displacements and 
compliance of DCB specimens reinforced with through thickness reinforcement known 
as Z-fibers. The Z-fibers were represented using discrete springs and a bi-linear cohesive 
damage law was used to represent Z-fiber failure as they are pulled out from a DCB 
specimen during testing. The closed-form solutions were used to calculate the rate of 
change of specimen compliance with delamination length. Linear elastic fracture 
mechanics was then used to calculate the strain energy release rate. Delamination growth 
was assumed to take place when the strain energy release rate became equal to the critical 
strain energy release rate of the parent laminate material. The analysis was coded in 
FORTRAN 90, yielding a dedicated analysis tool. Initial verification of the analysis was 
undertaken by comparing data from the analysis with experimental data from the 
literature. Agreement between analysis and experiment was reasonable. Data from the 
analysis was also compared with that of a previous analysis for predicting delamination 
growth in Z-fiber reinforced DCB specimens. Agreement between the current and 
original analysis methods was good. 
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Table 1. Details [17] of DCB specimen modeled in convergence study 


Arm height [mm], h 

1.61 

Specimen width [mm], b 

20.05 

Insert length [mm], a D 

49.6 

Flexural modulus [N/mm 2 ], E 

125500 


Table 2. Input data used in analysis of IM7/8552 based DCB specimen [12]. 


Laminate Details 

Z-fiber Details 

Bridging Law 

Ann height, h 

3 mm 

Diameter 

0.28mm 

Zo 

0.007mm 

Width, b 

1 8mm 

Length 

4mm 

Zf 

2mm 

Insert length 

73 mm 

Axial modulus 

227GPa 

k 0 

31446N/mm 

e d =e zd 

160GPa 

Spacing 

2mm 


G, c 

110J/m 2 

Areal density 

1.5% 


Table 3. Input data used in analysis of IMS/924 based DCB specimen with 0.5% areal 

density of Z-fibers [17]. 


Laminate Details 

Z-fiber Details 

Bridging Law 

Arm height, h 

1.615mm 

Diameter 

0.28mm 

Zo 

4.4 x 10" 3 mm 

Width, b 

20mm 

Length 

3.23 mm 

Zf 

1.615mm 

Insert length, a D 

50.7 mm 

Axial modulus 

227GPa 

k 0 

21638N/mm 

Arm length to 
Z-fiber region 

53.7 mm 

Spacing 

3.5mm 


E p =E zp 

126GPa 

Z-fibers / row 

5 

Gic 

268J/nf 

No of rows 

8 



Areal density 

0.5% 


Table 4. Input data used in analysis of IMS/924 based DCB specimen with 2% areal 

density of Z-fibers [17]. 


Laminate Details 

Z-fiber Details 

Bridging Law 

Arm height, h 

1 .63mm 

Diameter 

0.28mm 

Z 0 

4.4 x 10" 3 mm 

Width, b 

20mm 

Length 

3.26mm 

Zf 

1 ,63mm 

Insert length, a 0 

50.2mm 

Axial modulus 

227GPa 

k 0 

47165N/mm 

Arm length to 
Z-fiber region, 

a p 

51.6mm 

Spacing 

1.8mm 


E p =E zp 

126GPa 

Z-fibers / row 

11 

G, c 

268J/nr 

No of rows 

16 



Areal density 

2% 
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Figure 1. (a) Z-fibers held in foam preform positioned on laminate, (b) Ultrasonic 
hammer drives Z-fibers into laminate, (c) Completion of Z-fiber insertion. Preform 

material removed during insertion. 



Z-fibers misaligned relative to intended insert direction by 
an angle, 0. 


Figure 2. Cross-section of a laminate containing Z-fibers. 
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Figure 4. 2-Dimensional representation of Z-fiber reinforced DCB specimen. 


Specimen arm. 


(a) 


El 


F I i '' 

J^zp 1 zp \ ^ 


ki< 


he 


a 0 


aj- 


(b) 


Ri 

r 

© 

^ 

i 

X ■* 



1 


) 


Mn 


Figure 5. (a) Model of DCB specimen arm bridged by one Z-fiber row. (b) Free body 

diagram of the specimen arm. 
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Figure 6. Models of DCB specimen arms bridged by: (a) Two Z-fiber rows, (b) Three Z- 
libcr rows, (c) Arbitrary number of Z-fiber rows, n. 
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Figure 7. Idealized damage history of Z-fiber pull-out from a laminate. Illustration of 

single Z-fiber pull out test configuration. 
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Figure 8. (a) Schematic of DCB specimen with two different delamination lengths, 
(b) Load displacement response of DCB specimen with different delamination lengths. 
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Figure 9. Flow chart of solution procedure for predicting delamination growth in Z-libcr 

reinforced DCB specimens. 
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Figure 10. Load and displacement values as a function of dL/L, normalized by the 
corresponding experimental value, (a) Delamination length = 50mm. (b) Delamination 

length =99mm [17]. 


26 





Figure 1 1 . Specimen configuration used in Z-fiber reinforced DCB specimen 
investigation [12] (not drawn to scale). 



Figure 12. Estimated load displacement response of IM7/8552 DCB specimen containing 

1.5% areal density of Z-fibers. 
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Figure 13. Compliance calibration plot from analysis of IM7/8552 DCB specimen 
containing 1.5% areal density of Z-fibers. 



Delamination Length [mm] 

Figure 14. Modified beam theory plot from analysis of IM7/8552 DCB specimen 
containing 1.5% areal density of Z-fibers. 
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Figure 15. Estimated R-curve from analysis of IM7/8552 DCB specimen containing 


1 .5% areal density of Z-fibers. Comparison with ref 12. 



Crosshead Displacement [mm] 

Figure 16. Estimated load displacement response of IMS/924 specimen containing 0.5% 
areal density of Z-fibers. Comparison with ref. 17. 
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Figure 17. Estimated load displacement response of IMS/924 specimen containing 2% 
areal density of Z-fibers. Comparison with ref. 17. 
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Figure 18. Displacement vs delamination length of DCB specimen with 0.5% areal 
density Z-fibers. Comparison with ref. 17. 
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Figure 19. Displacement vs delamination length of DCB specimen with 2% areal density 

Z-fibers. Comparison with ref. 17. 


log C 

[log (mm/N)] 



-0 0.5% areal density 

-©- 2% areal density 

^developing 

2.36 

^developing 

1.96 

flsteady-state 

3.07 

^steady-state 

2.95 


- 0.5 . 


- 0.6 . 
- 0.7 . 


- 0.8 : 


- 0.9 : 


-1 . 

-1 1 Lj\ ; i ...j i I L... i 

0 N 1.70 1.75 1.80 1.85 1.90 


log a [log (mm)] 


1.95 



Figure 20. Compliance calibration plot from analysis of IMS/924 DCB specimen 
containing 0.5% and 2% areal density of Z-fibers. 


31 




J/3 


[mm/N] 


1/3 


0.8 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 


°20 



-H- 0.5% areal density 

-©- 2% areal density 

^Developing 

18.64-mm 

^Developing 

35.9-mm 

^Steady-state 

-1 .11 -mm 

^Steady-state 

2.68-mm 


100 


0 20 40 60 80 

Delamination Length [mm] 

Figure 2 1 . Modified beam theory plot from analysis of IMS/924 DCB specimen 
containing 0.5% and 2% areal density of Z-fibers. 
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Figure 22. Compliance calibration plot using test data from ref. 17 from tests on DCB 
specimens containing 0.5% and 2% areal density of Z-fibers. 
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Figure 23. Estimated R-curve from analysis of IMS/924 DCB specimen with 0.5% areal 
density of Z-fibcrs. Comparison with experiment. 



Figure 24. Estimated R-curve from analysis of IMS/924 DCB specimen with 2% areal 
density of Z-fibcrs. Comparison with experiment. 
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Figure 25. (a) DCB specimen modeled in Z-fiber row spacing study, (b) Bridging law for 
Z- fiber spacing case of 1.5-mm. (c) Bridging law for Z-fiber spacing case of 5-mm. 
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Figure 26. Load displacement response from analyses of DCB specimens with different 

Z-fiber spacing. 
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Figure 27. Load-point displacement versus delamination length curves from analyses of 
DCB specimens with different Z-fiber spacing. 
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